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The total hadronic cross-section cr 77 for the interaction of real photons, 77 — > hadrons, is extracted from 
a measurement of the cross-section of the process e + e~ — > e + e~7*7* — > (e + e~ + hadrons) using a luminosity 
function for the photon flux and form factors for extrapolating to Q 2 — 0. The data was taken with the OPAL 
detector at LEP at e + e~ centre-of-mass energies \/s ee ~ 161 GeV and 172 GeV. In the energy range 10 < W < 
110 GeV the total hadronic 77 cross-section <r 77 is consistent with the Regge behaviour of the total cross-section 
observed in 7p and hadron-hadron interactions. 



1. Introduction 

At high 77 centre-of-mass energies W = ^/s. 



77 



the total cross-section for the production of 
hadrons in the interaction of two real photons is 
expected to be dominated by interactions where 
the photon has fluctuated into an hadronic state. 



2. Kinematics 

The kinematics of the process e + c~ — > (e + e~ + 
hadrons) at a given \/s ee can described by the 
negative square of the four-momentum transfers, 
Qf — —qf, carried by the two (i = 1,2) incom- 
ing photons and by the square of the invariant 



Measuring the y/s dependence of the total mass of the hadronic final state, W 2 



hadronic 77 cross-section cr 77 should therefore 
improve our understanding of the hadronic na- 
ture of the photon and the universal high energy 
behaviour of total hadronic cross-sections. 

Before LEP the total hadronic 77 cross-section 
has only been measured for 77 centre-of-mass 
energies W below 10 GeV by PLUTO §, 
TPC/27 §] and the MD1 experiment g where 
the high energy rise of the total cross-section 
could not have been observed. Using LEP data 
taken at e + e~ centre-of-mass energies \/~s ee = 
130 - 161 GeV, L3 §| has demonstrated that 
the total hadronic 77 cross-section in the range 
5 < W < 75 GeV is consistent with the univer- 
sal Regge behaviour of total cross-sections. We 
present a measurement of the total hadronic 77 
cross-section in the range 10 < W < 110 GeV 
using OPAL data taken at v^cc = 161 GeV and 
172 GeV. 
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(?i + 92) 2 B- Events with detected scattered 
electrons (single-tagged or double-tagged events) 
are excluded from the analysis. This anti-tagging 
condition defines an effective upper limit Qm ax on 
the values of Q\ for both photons. This condition 
is met when the scattering angle 9' of the electron 
is less than the angle # max = 32 mrad between the 
beam axis and the inner edge of the acceptance 
of the detector or if the energy of the scattered 
electron is smaller than the minimum energy of 
35 GeV required for the tagged electron. 

3. Monte Carlo simulation 

The leading order (LO) QCD Monte Carlo gen- 
erators PYTHIA 5.722 || and PHOJET 1.05c § 
are used to simulate photon-photon interactions. 
PYTHIA is based on a model by Schuler and 
Sjostrand || and PHOJET has been developed by 
Engel based on the Dual Parton model (DPM) . 
The SaS-lD parametrisation of the parton distri- 
bution functions |10| is used in PYTHIA and the 
leading order GRV parametrisation |ll]] in PHO- 
JET. The fragmentation and decay of the parton 
final state is handled in both generators by the 
routines of JETSET 7.408 101. 
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4. Event selection 

Two-photon events are selected by requiring 
that the visible invariant hadronic mass, Wecal, 
measured in the electromagnetic calorimeter 
(ECAL), has to be greater than 3 GeV. At least 
3 tracks must have been found in the event and 
the sum of all energy deposits in the ECAL and 
the hadronic calorimeter (HCAL) has to be less 
than 45 GeV. The missing transverse energy of 
the event measured in the ECAL and the forward 
calorimeters (FD) has to be less than 5 GeV. No 
track in the event has a momentum greater than 
30 GeV/c. Finally, to remove events with scat- 
tered electrons in the FD or the silicon tungsten 
calorimeter (SW), the energy measured in the FD 
has to be less than 50 GeV and the energy mea- 
sured in the SW less than 35 GeV (anti-tagging 
condition). 

Additional cuts are applied to reject beam-gas 
and beam-wall background. On average the trig- 
ger efficiency for the lowest W range, 10 < W < 
30 GeV, is greater than 97% and it approaches 
100% for larger values of W. 

We use data corresponding to an integrated lu- 
minosity of 9.9 pb _1 at y/s cc = 161 GeV and 
10.0 pb^ 1 at Vi oc = 172 GeV. After applying all 
preselection cuts 55169 events remain. From the 
Monte Carlo (MC) it is estimated that after all 
cuts about 4% of all remaining events are e-f pro- 
cesses with Q 2 > 1 GeV 2 . The background from 
other processes apart from beam-gas and beam- 
wall interactions amounts to less than 1%. 



5. W reconstruction 

For measuring the total hadronic 77 cross- 
section cr 77 the value of W must be reconstructed 
from the hadronic final state. After the event 
selection a matching algorithm is applied in or- 
der to avoid double counting of particle momenta. 
The matching algorithm uses all the information 
of the ECAL, the HCAL, the FD and the SW 
calorimeters, as well as the tracking system. The 
four-momenta of the detected particles are used 
to calculate the visible invariant mass Wvis- 

A cut Wvi S > 6 GeV is applied to all prese- 
lected events. The W v i s distribution dN/dW V [ S 



measured at 172 GeV is shown in Fig. [l] where N 
is the number of selected events. They are well 
described by the MC simulations which have been 
normalized to the number of data events. 
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Figure 1. W V i S distribution for all selected events 
with Wvis > 6 GeV after background (bg) sub- 
traction at v^oc = 172 GeV compared to MC 
predictions. Statistical errors only are shown. 



6. Unfolding of the hadronic cross section 

The differential cross-section dcr ee /dW for the 
process e + e~ — ► (e + e~ -I- hadrons) has to be ob- 
tained from the Wvis distribution. The correla- 
tion between W v ; s and the generated invariant 
mass W for all selected PHOJET and PYTHIA 
events is shown in Fig. |^. The correlation is 
not very good due to hadrons which are emit- 
ted at small polar angles 8. These hadrons are 
either lost in the beam pipe or they are only 
detected with low efficiency in the electromag- 
netic calorimeters in the forward regions (FD and 
SW). The acceptance for PYTHIA is about 15% 
lower than for PHOJET at W = 40 GeV and it 
approaches the PHOJET acceptance of 60-65% 
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Figure 2. Correlation between W v i s and the gen- 
erated invariant mass W at y / s oc = 172 GeV for 
MC events. The vertical bars show the standard 
deviation (spread) of the Wvis distribution in each 
bin. For illustration purposes the points have 
been shifted by ±2 GeV. 
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Figure 3. The ratio of the number of selected 
events, A^ e i, to the number of generated events, 
iVg en , at a given generated invariant mass W at 
i/s ee = 172 GeV for MC events. The lower curves 
give this ratio for the diffractive and elastic events 
separately. 



for W > 80 GeV (Fig. g). The unfolding of 
these resolution effects, the correction for the de- 
tector acceptance and the background subtrac- 
tion is done by applying the unfolding program 
RUN p^| . The subtracted background does not 
include the remaining beam-gas and beam-wall 
interactions. Since the chosen bin size is not much 
larger than the resolution, bin-to-bin correlations 
are still sizeable. 

The differential cross-section dcr cc of the pro- 



hadrons) can be trans- 



lated into the cross-section cr 77 for the process 
77 — > hadrons using the luminosity function L 77 
for the photon flux § ,@ , @ 



der n 



dy 1 dQidy 2 dQi, 



a^(W,QlQ 2 2 ] 



d 4 L. 



77 



dy 1 dQjdy 2 dQl ' 



where y\ and yi denote the fraction of the beam 
energy carried by photons with y\yi ~ W 2 /s ee 
(neglecting Q 2 ). The cross-section for real pho- 
tons (Q 2 = 0) is derived by using appropriate 
form factors F(Q 2 ) which describe the Q 2 depen- 



dence of the hadronic cross-section: 

a 77 (W, Q 2 , Q 2 ) = F{Q\)F{Ql)a 11 {W, 0, 0) 

The luminosity function L 77 and the form factors 
F(Q 2 ) for the various W bins are obtained by 
applying the program PHOLUM 0. PHOLUM 
takes into account both transverse and longitudi- 
nally polarized photons. The uncertainty of the 
extrapolation to Q 2 = is estimated to be 5-7.5% 
by comparing the GVDM model with a simple p° 
form factor 0. This uncertainty is not included 
in the systematic error of the measurement. 

7. Systematic errors 

The two data samples at y/s cc = 161 GeV and 
172 GeV were independently analysed and the 
results for the total hadronic two-photon cross- 
section <t 77 are found to be in good agreement 
and are therefore averaged. 

Several distributions of the data are compared 
to the PYTHIA and PHOJET simulations af- 
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ter detector simulation in order to study whether 
the general description of the data by the MC 
is sufficiently good to use it for the unfolding of 
the cross-section. The MC distributions are all 
normalized to the data luminosity and the back- 
ground including the events with Q 2 > 1 GeV 2 
is subtracted from the data. 

In both MC models about 20% of the cross- 
section is due to diffractive and elastic events 
(e.g. 77 — * pp). This fraction is almost inde- 
pendent of W for W > 10 GeV. The selection 
efficiency for the diffractive and elastic events is 
very small and, although the rate is almost the 
same in both models, the selection efficiencies are 
very different. For a generated W of 50 GeV only 
about 6% of all generated diffractive and elastic 
events are selected in PYTHIA, whereas about 
20% are selected in PHOJET (Fig. §). Due to 
the small acceptance the detector correction has 
to rely heavily on the MC simulation for this class 
of events. 

Significant discrepancies are found in the distri- 
bution of the charged multiplicity n c h measured 
in the tracking chambers (Fig. ||). Both MC mod- 
els significantly underestimate the fraction of low- 
multiplicity events with n c h < 6 and overestimate 
the fraction of high- multiplicity events in compar- 
ison to the data. 

The energy -Efd measured in the forward de- 
tectors (FD) is shown in Fig. || for all selected 
events with Ep£, > 2 GeV. The good agreement 
of data and MC at large £pD shows that the back- 
ground from multihadronic Z° events and deep- 
inelastic e7 events is small and that this remain- 
ing background is reasonably well described by 
the MC 

Finally we plot the ratio P^/E v i s of the lon- 
gitudinal momentum vector Pl to the visible to- 
tal energy E vis (Fig- ||)- The ratio Ph/E V i S is 
peaked around 0.9 due to the Lorentz boost of 
the hadronic system. Data and MC are in rea- 
sonable agreement, different from the observa- 
tion in Ref . j| . Studies of beam-gas and beam- 
wall events show that most of these events have 
Ph/E vis > 0.85. The excess of the data over the 
MC seen at large P^/E v [ s is therefore consistent 
with 2% remaining background from beam-gas 
beam-wall events. 
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Figure 4. n c h distribution for all selected events 
with W V i S > 6 GeV after background subtraction 
at v^cc = 172 GeV compared to MC predictions. 
Statistical errors only are shown. 



Based on these observations the following sys- 
tematic errors are taken into account in the mea- 
surement of the total cross-section: 

• Both MC models describe the data equally 
well. We therefore average the results of 
the unfolding. The difference between this 
cross-section and the results obtained by us- 
ing PYTHIA and PHOJET alone are taken 
as systematic error. 

• For data and MC the cut on the charged 
multiplicity n c h was increased from n c h > 3 
to n c h > 5. This systematically shifts 
the cross-sections to lower values, mainly 
at small W . The unfolding was also re- 
peated using only data and MC events with 
3 < w c h < 9. This variation systematically 
increases the cross-section especially at high 
W where the average charged multiplicity is 
higher than at low W. The shifts are used 
as systematic errors. 

• The systematic error due to the uncertainty 
in the energy scale of the ECAL was esti- 
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Figure 5. The distribution of the energy in the 
forward calorimeter, -Efd, for all selected events 
with W V i S > 6 GeV after background subtraction 
at v / s cc = 172 GeV compared to MC predictions. 
Statistical errors only are shown. 
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Figure 6. The distribution of the ratio Pl/E v i s 
for all selected events with W V i S > 6 GeV after 
background subtraction at y/s cc = 172 GeV com- 
pared to MC predictions. Statistical errors only 
are shown. 



mated by varying the reconstructed ECAL 
energy in the MC by ±5%. 

• An overall normalisation uncertainty of 1% 
is due to the error on the luminosity mea- 
surement. 

• The lower limit on the trigger efficiency is 
taken into account by an additional system- 
atic error of 3% in the range 10 < W < 
38 GeV. 

• It is estimated that about 2% of the selected 
events could be due to beam-gas or beam- 
wall interactions. This value is therefore 
taken as additional systematic error. 

For the total error the statistical and the system- 
atic errors are added in quadrature. The values 
are given m Table |. 

8. Results 

The total hadronic cross-section a~ n for the 
process 77 — ► hadrons is shown in Fig. In 



the region W « 10 GeV the OPAL measurement 
is in agreement with the measurements at lower 
energies by PLUTO §, TPC/2 7 § within the 
experimental errors. 

The OPAL measurement shows the rise in the 
W range 10<IV<110 GeV which is characteris- 
tic for hadronic cross-sections. A similar rise was 
observed by the L3 experiment || , but the values 
for about 20% lower than the OPAL mea- 

surements. Several things should be noted which 
can explain parts of this discrepancy: First, the 
errors are strongly correlated between the W bins 
in both experiments. Secondly, L3 has used PHO- 
JET for the unfolding, whereas for the OPAL 
measurement the unfolding results of PHOJET 
and PYTHIA are averaged. The unfolded cross- 
section using PHOJET is about 5% lower than 
the central value. In both experiments the cross- 
sections obtained using PHOJET are lower than 
the cross-section obtained with PYTHIA. 

The total cross-section tr 77 is compared to sev- 
eral theoretical models. Based on the Donnachic- 
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W range (GeV) 


10- 16 


16 - 26 


26 - 38 


38 - 56 


56 - 80 


80 -110 


ct 77 [nb] 


385 


394 


398 


433 


485 


536 


stat. error 


± 4 


± 3 


± 4 


± 6 


± 8 


±15 


MC model 


±26 


± 4 


±20 


±21 


±21 


±43 


rich cut 


+ 5 


+ 13 


+ 24 


+ 45 


+ 74 


+ 95 


- 42 


- 37 


- 25 


- 21 


- 27 


- 34 


ECAL scale 


±37 


±28 


±13 


± 6 


± 6 


± 7 


luminosity 


± 4 


± 4 


± 4 


± 4 


± 5 


± 5 


trigger 


+12 


+12 


+ 12 








beam-gas 


- 8 


- 8 


- 8 


- 9 


-10 


-11 


total syst. 


+ 47 

- 62 


+ 34 

- 48 


+ 37 

- 36 


+ 5U 

- 32 


+ 78 

- 36 


+ 105 

- 56 


total error 


+ 47 

- 62 


+ 34 

- 48 


+ 37 

- 36 


+ 50 

- 32 


+ 78 

- 37 


+ 106 

- 58 



Table 1 

The total hadronic two-photon cross-section ct 77 and the contributions from the various systematic errors. 



Landshoff model |15]] , we test the assumption 
of a universal high energy behaviour of 77, 7p 
and pp cross-sections. The total cross-sections a 
for hadron-hadron and 7P collisions are well de- 
scribed by a Regge parametrisation of the form 

a = Xs e + Ys~ r ', (1) 

where ^fs is the centre-of-mass energy of the 
hadron-hadron or +p interaction. The first term 
in the equation is due to Pomeron exchange and 
the second term is due to Reggeon exchange Jl6) . 
The factors e = 0.0790 + 0.0011 and 77 = 0.4678 ± 
0.0059 are assumed to be universal and have been 
taken from Ref. jl6| together with the process de- 
pendent fit values of the parameters X and Y for 
the total hadronic +p and pp cross-sections. As- 
suming factorisation of the Pomeron term X, the 
total 77 cross-section can be related to the pp (or 
pp) and 7p total cross-sections at high centre-of- 
mass energies Vs 77 = V* 7P = v^pp where the 
Pomeron trajectory should dominate by 

cr 77 = — ^. (2) 
CTpp 

This simple ansatz gives a reasonable descrip- 
tion of <7 77 . Schuler and Sjostrand [^| give a to- 
tal cross-section for the sum of all possible event 
classes in their model of 77 scattering where the 
photon has a direct, an anomalous and a VMD 
component. They consider the spread between 
this prediction and the simple factorisation ansatz 



as conservative estimate of the theoretical band of 
uncertainty. We also plot the prediction of Engel 
and Ranft |7j which is implemented in PHOJET. 
It is in good agreement with the L3 measurement 
and significantly lower than the OPAL measure- 
ment. The steeper rise predicted by Engel and 
Ranft is in agreement with both measurements. 

9. Conclusions 

We have measured the total cross-section of the 
process 77 — * hadrons in the range 10 < W < 
110 GeV using the OPAL detector at LEP. 

Both MC models used fail to describe several 
distributions related to the hadronic final state 
like the charged multiplicity distribution. Further 
improvements of the description of the hadronic 
final state are necessary to reduce the systematic 
error of the measurement. It will also be impor- 
tant to gain a better understanding of the diffrac- 
tive and elastic processes for which the detection 
efficiency is found to be small. 

With the LEP2 data the high energy behaviour 
of the total 77 cross-section can be studied for the 
first time, extending the accessible W values by 
one order of magnitude up to W = 110 GeV. We 
observe the rise of the total 77 cross-section char- 
acteristic for the high energy behaviour of total 
hadronic cross-sections. A simple model based on 
Regge factorisation and a universal Donnachie- 
Landshoff fit to the total cross-sections of 7p and 
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Figure 7. function of W. The 

OPAL measurement is compared to the measure- 
ments by PLUTO §, TPC/2 7 [|, MD1 § and 
L3 Q. The inner error bars, which are some- 
times smaller than the symbol size, give the sta- 
tistical errors and the outer error bars the total 
errors. The data are compared to model predic- 
tions based on a Donnachie-Landshoff fit to total 
cross-sections JT5[ . The solid line gives the pre- 
diction using equation |^. The dotted line is the 
model of Schuler and Sjostrand ||. The model of 
Engel and Ranft used in PHOJET is shown as 
dashed line. 
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